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Abstract 

Prognostics, which deals with predicting remaining useful 
life of components, subsystems, and systems, is a key tech- 
nology for systems health management that leads to improved 
safety and reliability with reduced costs. The prognostics 
problem is often approached from a component-centric view. 
However, in most cases, it is not specifically component life- 
times that are important, but, rather, the lifetimes of the sys- 
tems in which these components reside. The system-level 
prognostics problem can be quite difficult due to the increased 
scale and scope of the prognostics problem and the rela- 
tive lack of scalability and efficiency of typical prognostics 
approaches. In order to address these issues, we develop 
a distributed solution to the system-level prognostics prob- 
lem, based on the concept of structural model decomposi- 
tion. The system model is decomposed into independent 
submodels. Independent local prognostics subproblems are 
then formed based on these local submodels, resulting in a 
scalable, efficient, and flexible distributed approach to the 
system-level prognostics problem. We provide a formulation 
of the system-level prognostics problem and demonstrate the 
approach on a four-wheeled rover simulation testbed. The re- 
sults show that the system-level prognostics problem can be 
accurately and efficiently solved in a distributed fashion. 

1. Introduction 

Prognostics is the process of predicting the end of (useful) life 
(EOL) and/or the remaining useful life (RUL) of components, 
subsystems, or systems. The prognostics problem itself can 
be divided into two distinct problems: (/) the estimation prob- 
lem, which determines the current state of the system, and (//) 
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the prediction problem, which, using the current system state 
estimate, computes EOL and/or RUL. In this paper, we focus 
on a model-based prognostics approach (Orchard & Vachtse- 
vanos, 2009; Daigle & Goebel, 201 lb; Saha & Goebel, 2009; 
Luo et al., 2008). In model-based prognostics, an underly- 
ing model of the system, its components, and how they fail 
is leveraged, where health state estimation is formulated as 
a joint state-parameter estimation problem, typically using a 
filtering approach, and prediction is formulated as a simula- 
tion problem (Daigle, Saha, & Goebel, 2012). 

To the best of our knowledge, all prognostics research to date 
has been focused on individual components, and determining 
their EOL and RUL, e.g., (Orchard & Vachtsevanos, 2009; 
Saha & Goebel, 2009; Daigle & Goebel, 2011a; Celaya et 
al., 2011; Bolander et al., 2010; Luo et al., 2008; Bying- 
ton et al., 2004). However, in many cases, the desired infor- 
mation is the EOL of the system, which is obtained through 
system-level prognostics. Generally, the EOL of a system de- 
pends on its constituent components and how they interact. 
Approaching this problem from the centralized perspective 
becomes very difficult, as common (centralized) prognostics 
algorithms may not scale to the system level. 

In order to address the problems with centralized approaches, 
in recent work, we have developed a distributed model-based 
prognostics architecture that allows the decomposition of a 
large prognostics problem into several independent local sub- 
problems from which local results can be merged into a global 
result (Daigle et al., 2011; Daigle, Bregon, & Roychoudhury, 
2012). Since each local subproblem can be solved indepen- 
dently, each can be assigned to a different processing unit 
and be solved in parallel. Such a distributed approach is in 
contrast to other proposed distributed prognostics architec- 
tures in which the global problem is not decomposed and 
the computation is distributed onto multiple processing units, 
e.g., (Saha, Saha, & Goebel, 2009). Our distributed approach 
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scales well and the resulting subproblems are typically small 
and easy to solve, resulting in an efficient and flexible dis- 
tributed solution to the prognostics problem. Such an ap- 
proach has obvious advantages when applied to the system- 
level prognostics problem. In this paper, we formulate the 
system-level prognostics problem and propose a solution us- 
ing this distributed prognostics framework. We apply our 
system-level prognostics approach to a rover testbed and pro- 
vide results in simulation to empirically demonstrate and val- 
idate the approach. 

The paper is organized as follows. Section 2 formulates 
the system-level prognostics problem and overviews the pro- 
posed distributed solution. Section 3 describes the estima- 
tion problem, and Section 4 describes the prediction problem. 
Section 5 presents the rover case study, and shows prognos- 
tics results in simulation. Section 6 concludes the paper. 

2. System-level Prognostics 

While most prognostics approaches focus on individual com- 
ponents, in most practical cases it is actually the EOL of the 
system that must be determined. With this prediction, the 
future usage of the system may be optimally planned to max- 
imize system life and to schedule system-wide maintenance 
activities. It is often important to take a system-level per- 
spective of prognostics, because the degradation of individual 
components is often coupled, i.e., the way one component de- 
grades is dependent on how a connected component degrades. 
This may occur, for example, if one component provides the 
inputs to another component, in which case, prognostics of 
the latter component cannot be performed in isolation. 

In this section, we first define the system-level prognostics 
problem. We then introduce the system-level prognostics ap- 
proach and architecture using a distributed prognostics frame- 
work that is based on structural model decomposition. 

2.1. Problem Formulation 

The goal of system-level prognostics is the prediction of the 
EOL and/or RUL of a system. We assume the system model 
may be generally defined as 

x(t) = f(t,x(t),0(t),u(t),v(t)), 
y (t) = h(t,x(t),0(t),u(t),n(t)), 

where x(t) E R nx is the state vector, 0(t) E R n<? is the 
unknown parameter vector, u (t) E R nu is the input vector, 
v(f.) E R nv is the process noise vector, f is the state equation, 
y (t) E R ny is the output vector, n(i t) E R n " is the measure- 
ment noise vector, and h is the output equation. 1 This model 
describes both the nominal behavior and faulty behavior, in- 
cluding the fault progression functions. 


’Here, we use bold typeface to denote vectors, and use n a to denote the 
length of a vector a. 


In system-level prognostics, we are interested in when the 
performance of a system lies outside some desired region 
of acceptable behavior. The desired performance is ex- 
pressed through a set of n c constraints, Ceol = 
where c* : R Ux x R n * x R nu — » B maps a given point 
in the joint state-parameter space given the current inputs, 
(x(£), 9(t). u (£)), to the Boolean domain B = [0, 1], where 
Ci(x(£),0(£),u(£)) = 1 if the constraint is satisfied. If 
d(x(t), 9(t), u(t)) = 0, then the constraint is not satis- 
fied, and the behavior of the system is deemed to be un- 
acceptable. These deterministic constraints may refer to 
component-level, subsystem-level, or system-level specifica- 
tions or requirements and define a fixed partition of the state- 
parameter-input space into acceptable and unacceptable re- 
gions of behavior. When the constraints are violated, it does 
not necessarily refer to a hard failure, but any point at which 
the operational risk is too large to continue system operation, 
or future behaviors of the system will be in some way unac- 
ceptable. At this point we say the system has no useful life 
remaining. 

These individual constraints may be combined into a single 
system-level threshold function Teol ' R Ux x R n * x R Uw — > 
B, defined as 

T E OL(x(t),0(t),u(t)) = 

1. 0 € {Ci(x(t), 0(0,11^))}^ 

0, otherwise. 

Teol evaluates to 1, i.e., the system has reached an unac- 
ceptable region of behavior, when any of the constraints are 
violated. EOL is then defined as 

EOL(t P ) = 

inf {t € R : t > t P A T£ol(x(<), 0{t), u (t)) — 1}, 

i.e., EOL is the earliest time point at which Teol is met (eval- 
uates to 1). RUL is expressed using EOL as 

RUL(t P ) = EOL(tp) -t P . 

Note that because x(t) is a random variable, EOL and RUL 
must necessarily be random variables also. 

2.2. Prognostics Approach 

In order to make an EOL or RUL prediction for the system, 
the initial state from which to make a prediction is required. 
In general, this initial state is not directly observed, and must 
be estimated. Therefore, there are two sequential problems 
for prognostics: the estimation problem and the prediction 
problem. The estimation problem is to find a joint state- 
parameter estimate p(x(t),0(t) |y 0:t ) based on the history of 
observations up to time £, yo : *. This estimate is represented 
as a probability distribution because, generally, the system 
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state is not directly observed, and there is sensor noise, n(£), 
and process noise, v(£). At a given prediction time, tp, the 
prediction algorithm uses the joint state-parameter estimate 
p{x{tp),0(tp) |y 0 :(p) and computes p(EOL(tp)\y 0 -.t P ) and 
p(RUL(tp)\y 0 : t p ). Along with the uncertainty in the state- 
parameter estimate, process noise and uncertainty in the fu- 
ture inputs to the system all contribute to the uncertainty in 
the EOL/RUL prediction. 

This system-level prognostics problem, consisting of estimat- 
ing the system state and then predicting its evolution to EOL, 
can be solved using component-level approaches by treating 
the entire system as a single component and applying these 
approaches directly. However, for a large system, both the es- 
timation and prediction problems are correspondingly large. 
Due to the large state-parameter dimension, a centralized ap- 
proach does not scale well, and can be very inefficient. 

Therefore, we propose to decompose the global system- 
level prognostics problem into independent local subprob- 
lems, such that the solutions to the local subproblems may 
be easily merged to form the solution to the global prognos- 
tics problem. This forms a naturally distributed approach in 
which the local subproblems, since they are independent, may 
be solved in parallel, thus providing scalability and efficiency. 
Further, the approach allows different algorithms to be em- 
ployed on each subproblem. The subproblems often corre- 
spond directly to component-level prognostics problems, and 
the approach provides a mechanism to combine component- 
level prognostics results into system-level results. 

In (Daigle et al., 2011), we developed such a distributed solu- 
tion to the estimation part of the prognostics problem, based 
on the concept of structural model decomposition (Pulido & 
Alonso-Gonzalez, 2004). In recent work, the same concept 
was used to decompose the prediction problem (Daigle, 
Bregon, & Roychoudhury, 2012). Structural model decom- 
position allows one to decompose a system model into a set 
of submodels for which local prognostics problems can be 
directly defined. The global model of the system, denoted as 
M , is defined as follows. 

Definition 1 (Model). The model of a system, Al, is a tuple 
M = ( X, ©, U, y, C), where X is the set of state variables 
of x, O is the set of unknown parameters of 0 , U is the set of 
input variables of u, Y is the set of output variables of y, and 
C is the set of model constraints of f, h, and Ceol- 

Informally, a model consists of a set of variables and a set 
of constraints among the variables. While technically f and 
h themselves are (complex) constraints, we represent them 
instead as sets of simple constraints. This view is also more 
consistent with the way modelers describe f and h, i.e., as sets 
of equations, each describing a single state or output variable. 

Model decomposition is accomplished by assigning some 
variables as local inputs for which the values are known (e.g.. 


they are directly measured). In this way, the submodels are 
made computationally independent of each other. Within this 
scheme, a submodel is then defined as follows. 

Definition 2 (Submodel). A submodel Mi of a system model 
M = (X,G,U,Y,C) is a tuple Mi = 
where X { C X, ©< C ©, U { C X U U U Y, and Y { C Y are 
the state, parameter, input, and output variables, respectively, 
and Ci C C are the submodel constraints. 

For distributed prognostics, we find a set of submodels that 
satisfy a certain set of properties. For distributed estimation, 
the submodels use Ui C U U (Y - Yi), and we find a set of 
minimal submodels such that each Yi is a singleton, and over 
all Y i: Yj where i ^ j, Yi PI Yj = 0. So, each submodel uses 
some global model inputs and some measured values as lo- 
cal inputs, and, in this way, the submodels become decoupled 
and may be computed independently from each other. By cre- 
ating submodels with one output variable each, we maximize 
the number of estimation submodels and the opportunity for 
parallelization of the estimation task. By making the sub- 
models minimal, they require no constraints or variables that 
are not strictly necessary to compute Y{. An algorithm for 
computing the set of submodels with these properties is given 
in (Daigle et al., 2011), which is based on the model decom- 
position algorithms presented in (Pulido & Alonso-Gonzalez, 
2004; Bregon, Biswas, & Pulido, 2012). 

For distributed prediction, the submodels use Ui C Up , 
where Up C X U U. Here, Up is a set of variables whose 
future values can be hypothesized. In the centralized case, 
Up = U. We find a set of minimal submodels such that 
each submodel has at least one c G Ceol belonging to Ci , 
and over all submodels, Ceol is covered. This ensures that 
Teol ™ay be computed for the system; since Teol is 1 
whenever any of the constraints in Ceol are violated, we can 
independently evaluate when those individual constraints will 
be violated and then take the minimum to obtain the system 
EOL. An algorithm for computing the set of submodels with 
these properties is given in (Daigle, Bregon, & Roychoud- 
hury, 2012). Both decomposition algorithms work in a sim- 
ilar way; essentially, they start with a variable or constraint 
that must be computed in the local submodel, and then trace 
the dependencies backwards until local inputs are reached, 
including all variables and constraints found throughout the 
search within the submodel. 

Note that the problem of defining Up is critical to obtaining 
accurate results for system-level EOL in a distributed manner. 
On average, the most accurate result will be achieved when 
the system model is directly used for prediction, because it 
captures all the interdependencies between the components. 
In the general case, damage could be progressing in multiple 
components, and damage progression in one component may 
have an effect on damage progression in another component 
due to their coupling. In such cases, for system-level prog- 
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Prognostics 



Figure 1. Sample system-level prognostics architecture. 


nostics the components cannot be decoupled due to these in- 
teractions, and the prediction problem cannot be decomposed 
into two independent problems, one for each component. It is 
only appropriate to neglect these interactions when they are 
either negligible or predictable a priori. It will be shown in 
Section 5 how this is an important consideration. 

2.3. Prognostics Architecture 

A sample system-level prognostics architecture based on the 
distributed framework is shown in Fig. 1. In discrete time 
k, and using a discrete-time version of the model, the dam- 
age estimation module takes as input both u k and y* and 
splits them into local inputs and outputs for the submodels. 
Estimation is performed for each submodel using an appro- 
priate algorithm, computing local state-parameter estimates 
p(x x k ,Q\ |yQ. fc ). Some of these local estimates are merged 
corresponding to the prediction submodels. For example, 
submodel builds its local state using the estimates from 
the estimators of Ad 2 , M 3 , and Ad 4 . The local predictors 
compute local EOL/RUL predictions p(EOL l kp \y l 0:kp ) and 
p(RUL l kp \y l 0:kp ) at given prediction time kp based on the 
local EOL constraints. Local predictions are then merged into 
global predictions p(EOL kp |y 0: * P ) and p(RU L kp \y 0:k p) 
by taking the minimum of the local predictions. 

3. Distributed Estimation 

As described in Section 2, in our distributed estimation 
scheme, the local estimator for each submodel Mi produces 
a local estimate p(x k . OHy^.^), where x k C x k , 0\ C 0 k , 
and y k C y k . Here, the local inputs used, u\ consist of ele- 
ments from both u and y, where measured values are directly 
used as local inputs. The estimation problem is decomposed 
by finding a set of minimal submodels that together cover the 
subset of x and 6 required for prediction, by using these local 
inputs. This approach to distributed estimation is different 
from approaches like the distributed decentralized extended 
Kalman filter (Mutambara, 1998) or other estimation fusion 


techniques (Sinha et al., 2008) where local estimates are com- 
municated between local estimators. Here, local estimators 
do not communicate and operate completely independently. 

In order to effectively perform joint state-parameter estima- 
tion, the system should be observable, among other require- 
ments. If the global model is structurally observable, then we 
are guaranteed that the local submodels for estimation are as 
well (Moya et al., 2010). 

Any suitable algorithm may be used for joint state-parameter 
estimation. In this paper, we use an unscented Kalman fil- 
ter (UKF) (Julier & Uhlmann, 1997, 2004) with a variance 
control algorithm (Daigle, Saha, & Goebel, 2012). The UKF 
assumes the general nonlinear form of the state and output 
equations described in Section 2, but restricted to additive 
Gaussian noise. 

We summarize the main details of the UKF below, and refer 
the reader to (Julier & Uhlmann, 1997, 2004) for details. In 
the UKF, distributions are approximated using the unscented 
transform (UT). The UT takes a random variable x £ R nx , 
with mean x and covariance P xx , that is related to a second 
random variable y £ R ny by some function y = g(x), and 
computes the mean y and covariance P yy using a minimal set 
of deterministically selected weighted samples, called sigma 
points (Julier & Uhlmann, 1997). X 1 denotes the zth sigma 
point from x and w l denotes its weight . 2 The sigma points 
are always chosen such that the mean and covariance match 
those of the original distribution, x and P xx . Each sigma 
point is passed through g to obtain new sigma points y, i.e., 

y = g(*‘) 


2 Sigma point weights do not directly represent discrete probabilities, so are 
not restricted to [0, 1]. 
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with mean and covariance calculated as 


step is: 


P 


y = 

i 

«v=Y, wi (y i -y)(y i -y) T - 

i 


In this paper, we use the symmetric unscented transform, in 
which 2 n x + 1 sigma points are symmetrically selected about 
the mean according to (Julier & Uhlmann, 2004): 


w 


(jlx 

, 2 (n x + At) 

( _ 
x, 


* = 0 

i — 1, . . . , *2iTix 

?: = o 


x l = 


x+ ^(ns+KjP • • • > ^X 

X (Tlx Iij } i = “h 1 f • • • j 271 j; , 

where ^>/(n x 4- k)P xx ^ refers to the ?'th column of the ma- 


trix square root of ( n x -f Av)P xx . Here, k is a free parameter 
that can be used to tune higher order moments of the distribu- 
tion. If x is assumed Gaussian, then selecting n = 3 — n x is 
recommended (Julier & Uhlmann, 1997). 


In the filter, first, n s sigma points are derived 

from the current mean and covariance estimates 

Pk-i\k-i using a sigma point selection algorithm. The pre- 
diction step is: 


~ X * % 

Xk\k-\ — = U . . . ,n s 

y'k\k~i = = i .•••.«* 

n a 

*‘k\k - 1 = yyx k \^ 

i 

n a 

Yk\k-1 = ^2 Wt ^ k \k-l 


with 


Pk\k-1 = Q+ 

n a 

-*k\k-l)(Xi\k-l-*k\k-l) T , 

i 

where Q is the process noise covariance matrix. The update 


n s 

P yy =K 4- - yfc|A;-l)(3 ; fc|Ar~l ~ Yk\k-l) T 

i 

n a 

p *v = Yl wi ^ x k\k-i — *fcifc-i)(y3b|jb-i -yM*-i) T 

i 

K fc = P^P- 1 

Xfc|A: = Xk\k-1 + Kfc(yfc ~ Yk\k-l) 

Pk\k = Pk\k-1 - KfcPyyKj^, 

where R is the sensor noise covariance matrix. 

Joint state-parameter estimation is accomplished in the UKF 
by augmenting the state vector with the unknown parameters, 
and the corresponding diagonal elements of the process noise 
matrix, Q, are set to nonzero values. In this way, the param- 
eter estimates become time-varying and are modified by the 
filter using the measured outputs. 

The variance values in Q associated with the unknown pa- 
rameters determine both the rate of parameter estimation con- 
vergence and the estimation performance once convergence 
is achieved, therefore, techniques have been developed to 
tune this value online to maximize performance, e.g., (Liu & 
West, 2001; Orchard, Tobar, & Vachtsevanos, 2009; Daigle, 
Saha, & Goebel, 2012). We adopt the approach presented 
in (Daigle, Saha, & Goebel, 2012), in which the algorithm 
tries to control the variance of the hidden wear parameter 
estimate to a user-specified range by modifying the process 
noise variance. Effectively, the algorithm increases the vari- 
ance when the relative parameter spread is below the desired 
level, and decreases it otherwise. With the proper settings, 
the parameter estimates converge quickly and track with high 
accuracy and precision. 

4. Distributed Prediction 

Each local prediction module takes as input local state- 
parameter estimates formed from the local estimators, as dis- 
cussed in Section 2. The required estimates must be con- 
structed from the local estimates of the submodels used for 
estimation. A prediction submodel has a set of states Xi 
and parameters 0 j, and it must construct a local distribu- 
tion p(x\,O l k |yo : A:)* To do this, we assume that the local 
state-parameter estimates may be sufficiently represented by 
a mean fj, 1 and covariance matrix £*. For each prediction 
submodel Mu we combine the estimates from estimation 
submodels that estimate states and parameters in Xi U 0* 
into /. l 1 and covariance E*. If there is overlap in the state- 
parameter estimates, i.e., if two submodels both estimate the 
same state variable x or parameter 0 , then this may be re- 
solved by a number of techniques, e.g., taking the estimate 
with the smallest variance, or taking an average. Note that, 
due to the decomposition into independent local submodels. 
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Algorithm 1 EOL Prediction 


Inputs: {(x^ ) ,e^ ) ),u/^ ) }f =1 

Outputs: [EOL^.u^yfLi 

for j = 1 to N do 

k 4— kp 


4 (j) 

9 


*(j) 

x k P 

e kp 


Predict uj. 

while T' EOL {x k (j \ uj.) = 0 do 


Predict u*. 

p{0\ 


Q‘U) 

"fc+l 

iU) 
H + 1 




-pWUiK 


Hi) a Hi) 


k 4 — k T 1 
iti) , i(J ) 

^k + \ 

e\ U) <- ei U) 

end while 




'k+1 


EOL 

end for 


* 0 ) 

k P 


we recover only an approximation to the joint posterior dis- 
tribution as would have been found by a global estimator. In 
particular, covariance information is lost due to the decou- 
pling and will appear as zeros in the merged covariance ma- 
trix. As shown in (Daigle et al., 201 1) and as will be seen in 
Section 5, the approximation still results in accurate predic- 
tions. 

Given the mean and covariance information, we represent the 
distribution with a set of sigma points derived using the un- 
scented transform. Then, as in (Daigle & Goebel, 2010), each 
sigma point is simulated forward to EOL, and we recover the 
statistics of the EOL distribution given by the sigma points. 

The prediction algorithm is executed for each submodel i, de- 
riving local EOL predictions using its local threshold func- 
tion based on the local EOL constraints. The pseudocode for 
the prediction procedure is given as Algorithm 1 (Daigle & 
Goebel, 201 lb). For a given submodel Mi , each sigma point 
j is propagated forward until Tfe OL (x l £*\ 0^ j) ) evaluates to 
1. The algorithm hypothesizes future inputs u*.. 

Each prediction submodel Mi computes a local EOL/RUL 
distribution, i.e., p{EOIS kl ,\ y’ 0:fcj> ) and p(RUL i kp |yj :fcp ). 
The system EOL is determined by the minimum of all the 
local distributions, since Teol of the system is 1 whenever 
any of the local constraints are violated, and each local distri- 
bution is associated with a subset of these constraints. Specif- 
ically, for m prediction submodels, 

p(EOL kp \y 0:kp ) = min({ P (£;OL* fcp |yj :fcp )}- 1 ). 

To compute this, we sample from each local EOL distribution 
and take the minimum of the local samples. This is repeated 
many times and the statistics of the global EOL distribution 
are computed. 


5. Case Study 

In this section, we apply our system-level prognostics ap- 
proach to a four-wheeled rover testbed developed at NASA 
Ames Research Center. We develop a model of the rover, and 
demonstrate the approach using simulated scenarios. 

5.1. Rover Modeling 

The rover model was originally presented in (Balaban et al., 
2011). In this section we summarize the main features and 
include some extensions to the model. 

The rover consists of a symmetric rigid frame with four 
independently-driven wheels. The wheel speeds are governed 
by 

UFL = "7 ( TmFL ~ TfFL ~^glFL + T^FI,) (Cl) 

Jfl 

UFR = ~r ( TmFR - TfFR ~ T giFR - T grFR ) (c 2 ) 

Jfr 

UBL = "7 ( TmBL ~ r fBL ~ r glBL + TgrBL) (C 3 ) 

Jbl 

UBR = ~r ( r mBR ~ r fBR ~ r glFR ~ r g rBR) (c 4 ) 

Jbr 

The F, B, L, and R subscripts stand for front , left , 
back, and right, respectively. Here, for wheel w £ 
{FL, FR, BL, BR }, J w denotes the wheel inertia; r mw = 
kri w is the motor torque, where i w is the motor current and 
kr is an energy transformation gain; rj w = fj,f w uj w is the 
wheel friction torque, where f.ij w is a friction coefficient; 
Tgiw = rwVgis (v w — v) is the torque due to slippage, where 
r w is the wheel radius, fi g i s is a friction coefficient, v w is the 
translational wheel velocity, and v is the translation velocity 
of the rover body; and r grw = r w fj, grw uj cos 7 is the torque 
due to the rotational movement of the rover body, where f. i grw 
is a friction coefficient, u is the rotational velocity of the rover 
body, and 7 = arctanZ /6 with / being the rover length and 6 
being its width. 

We consider here friction-based damage progression in the 
motors, resulting in an increase in motor friction over time, 
which will lead to an increase in power consumption. For 
wheel w , is governed by (Daigle & Goebel, 2011b) 


A/FL = VfFL HfFL U FL 

(C 5 ) 

f'fFR = Vjfr M/FB U RR 

(C 6 ) 

P'fBL = "/BL VfBL U BL 

(c 7 ) 

l 1 f BR = VfBR fJ'fBR U R r, 

(Cs) 


where for wheel w, v f w is an unknown wear coefficient. 

The translational velocity v of the rover is described by 

V = — (FglPL + Fgl FR + FglBL + F gi & R ) , (C9) 
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Figure 2. Rover forces. 


where m is the rover mass, and for wheel w y Fglw = 
Vgiwiv w — v) is the force due to slippage. The rotational 
velocity lj is described by 

u = j(dcos'yF glFR + dcos'yF giBR -d cos'? F g i FL 

— d cos 7 F g iBL — dF grFL — d.F grFR — dF grB L 

— dF g rBR)’ (Cio) 

Here, J is the rotational inertia of the rover, d is the distance 
from the center of the rover to each wheel, and for wheel w y 
F grw — HgrivU is the force due to the rotational movement of 
the rover body. The rover forces are illustrated in Fig. 2. 

The wheels are driven by DC motors with PI control that sets 
the voltages V applied to the motors. The motor currents are 
governed by 


The temperature of the windings for the motors are governed 
by 


TdFL = (^FlF - h dF i(TdFL - Tm F L)) 
TdFR = jr (i 2 frF — h d FR{T d FR — T mFR )) 

FdBL = Y d BL R - hdBL(T dBL - T m £ L )) 

TdBR = -J~ d (* 2 BrR — hdBR(TdBR - T mBR )) , 


(C19) 

(c 2 o) 

(C 2 l) 

(C22) 


where J d is the thermal inertia of the windings, and for wheel 
w y h dw is a heat transfer coefficient, and T mw is the motor 
surface temperature. It is assumed that heat is lost only to 
the motor surface, and that winding resistance R is approx- 
imately constant for the temperature range considered. The 
surface temperature of the motor for wheel w is given by 


TmFL = ~r(hdFL(T d FL ~ T m FL) - /l a FL(T m FL — T a )) 
Js 

(c 23) 

TmFR = -J-(hdFR{TdFR ~ T m FR) ~ h a FR{T mF R — T a )) 
J 8 

(c 24) 

T m BL = ~r{hdBL{TdBL ~ T m BL.) ~ h a BL(T m BL ~ T a )) 
Js 

(C25) 

TmBR ~ -r(hdBR(T d BR ~ TmBR.) ~ h aBR {T m BR ~ T a )), 
Js 

(^26) 


IFL = J^FL - i FLRFL - k u U FL ) (Ql) 

IFR = jr(VFR - IFRRFR “ k^UJpR) (^12) 

IBL = J^(YbL ~ iBLRBL - kuWBL.) (^ 13 ) 

IBR = VBR ~ IBrRbR ~ kujUJBR)- (Cl4> 


Here, L is the motor inductance, k u is an energy transforma- 
tion term, and for wheel w , R is the motor resistance. The 
voltages applied to the motors are determined by the con- 
trollers, where for wheel w y V w = P* (u w — c j w ) + I * e^, 
where P is a proportional gain, u w is the commanded wheel 
speed, / is an integral gain, and ei w is the integral error term. 


The integral error terms are governed by 


&iFL = U F l — U FL 

(CIS) 

eiFR — u FR — LOpR 

( c 16) 

CiBL ~ UBL ~ WBL 

(c 17) 

&iBR = UbR — WbR- 

(^18) 


The motor windings heat up as current passes through them. 


where J s is the thermal inertia of the motor surface, and for 
wheel w y h aw is a heat transfer coefficient, and T a is the am- 
bient temperature. Heat is transferred from the windings to 
the surface and lost to the environment. 

The batteries, which are connected in series, are described 
by a simple electrical circuit equivalent model that includes a 
large capacitance Cb in parallel with a resistance Rp y together 
in series with another resistance i? s . 3 The battery charge vari- 
ables qi are governed by 

<7i = —Vi/Rpi - (}fl + ifr 4- ibr 4- Ibl) (C 27 ) 

Q2 = —V 2 / Rp2 - (}FL + *FF + ^BR 4- Irl) (c 2 s) 

<73 = -Vz/Rpz - (}fl -I- ^fr 4 ibr 4- ibl) (Q 29 ) 

</4 = —V4 / Rp 4 - ( iFL + 'IFR + IBR + (C30) 


3 We use a simple model here only for demonstration purposes. More detailed 
battery models for prognostics can be found in the literature, e.g., (Saha & 
Goebel 2009). 
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Table 1. Estimation Submodels 


The available sensors measure the voltages of the batteries, given as 


V\ = Ql/Cbl - Rs 1 * (*FL + IFR + iBR + i BL ) foil) 

V 2 * = 92 / C&2 - Rs2 * (iFL + IFR + iBR + ^Bl) (C 32 ) 

^ 3 * = 93/^63 — #s3 * (*FL + ^FR + iBR + *Fl) (^ 33 ) 

Kf = q^/CbA - Rs4 * (*FL + *FR + IbR + ^Bl), (^34) 

the battery current, 

h — i FL H- i>FR + iBR + IBL, (C35) 


the motor currents, 

*FL = ? 'FL 
^FR = i FR 
i*BL ~ ^BL 
* BR = ^BRi 

the motor surface temperatures, 

T mFL = TmFL 
T mFR = TmFR 
TniBL ~ TmBL 
T mBR = TmBR, 

and the wheel speeds, 

u fl = wfl 
w FR = U FR 
U*BL = U BL 
U*BR = “BR- 


(c 3 e) 

(C37) 

(C38> 

(C39) 


(C40) 

(C41) 

(C42) 

(C43) 


(C44) 

(C45) 

(^46) 

(C 47 ) 


Vi > V~ (C48) 

V 2 > V (C49) 

V 3 > V (c 50 ) 

V 4 > V~, (c 5 i) 

where the voltage threshold is given by V~ = 9.6 V, and for 


battery?:, Vi = qi/Cbi-R S i*(i F L+iFR+iBR+iBL)- We are 
also interested in when the motor temperature gets too high. 
The motor windings are designed to withstand temperatures 
up to a certain point, after which, the insulation breaks down, 
the windings short, and the motor fails (Balaban et al., 2010). 
The constraints are given as 

TmFL < T + (C 52 ) 

T mF R < T + (C53) 

T m BL < T + (C54) 

TmBR < T +, (C55) 

where the temperature limit is given by T+ = 70° C. The 
rover cannot be operated when any of these constraints, c 4 $- 
C55, are violated. 

In the general case, we consider uncertainty in the friction 
wear parameters is/fl , v/fr, v/bl> an d v/br\ the heat 
transfer coefficients hdFL> hdFR> hdBL > ^fl» ^oFb» 
h a BL, and h a BR\ the battery capacitances Cb 1 , Cb 2 , Cb 3 , and 
CbA\ and the battery resistances R s 1 , R s2 , R s 3 , and R s4 . Sen- 
sor and process noise were estimated based on data from the 
actual rover testbed. 


Here, the " superscript indicates a measured value. 

We are interested in predicting when any of the rover batter- 
ies are at their voltage threshold, beyond which the batteries 
will be damaged (Saha & Goebel, 2009). The constraints are 


5.2. Results 

To demonstrate the validity of the approach, we describe 
two scenarios for system-level prognostics of the rover. In 
the first, the rover is operating nominally without any faults 
present, and in the second, friction damage is progressing 


8 


Annual Conference of the Prognostics and Health Management Society 2012 


Submodel 

~Xi 


" 77 T" 

~Y~ 

~CT~ 

M 17 

<71 

L’ti , Rsi 

iFL,iFR,lBL,iBR 

0 

C27,C48 

M 18 

Q2 

Cb2i Rs2 

IFL , IFR, IBL , IBR 

0 

C281C49 

M 19 

Q3 

C&3, Rs3 

iFL,iFR,lBL,lBR 

0 

C29,C50 

M20 

<74 

Cb4, Rs4 

lFL,iFR, lBL,iBR 

0 

C 30 .C 51 

M 21 

TdFL.TmFL 

hdFL,haFL 

IFL 

0 

Cl9,C23*C52 

M 22 

TdFR, TmFR 

hdFR, haFR 

IFR 

0 

C2O.C24.C53 

M 23 

TdBL,TmBL 

hdBL,haBL 

\BL 

0 

C2i.C25.C54 

M 24 

TdBRjTmBR 

hdBR . h a BR 

i>BR 

0 

C22.C26.C55 


Table 2. Prediction Submodels Using Motor Currents as Local Inputs 


on one motor. In both cases, the rover travels between var- 
ious waypoints, moving at an average speed of 0.5 m/s. The 
unknown parameters are initialized incorrectly (with around 
10% error) so that the local estimators must converge to the 
true values. In both cases, the estimation step is performed 
in a distributed manner using the set of submodels derived 
by using measured values as local inputs, shown in Table 1. 
For example, submodel Mi computes an estimate of V{ us- 
ing the measured value of i b as a local input, and using the 
minimal set of constraints to do this. For the prediction sub- 
models, as will be shown, the correct submodels to use de- 
pends on the scenario, and illustrates when and when not the 
prediction step can be decomposed. 

5.2.1. Nominal Operation 

We first consider a scenario involving nominal, fault-free op- 
erations. In this case, EOL occurs around 3 h. An RUL pre- 
diction is made every 500 s. With the rover traveling at an 
average speed of 0.5 m/s, the motor currents average to about 
0.15 A each and so the total current draining from the four 
batteries is 0.6 A. Since these values do not vary much dur- 
ing nominal operation, we can use the motor currents as local 
inputs for the model decomposition. These submodels are 
shown in Table 2. Note that the estimates from the estima- 
tion submodels M\-M$, are used directly in the prediction 
submodels M 17 -M 24 , respectively, and that estimation sub- 
models M$-M\q, are not necessary. Note also that the pre- 
diction submodels do not compute any outputs, rather, their 
goal is to compute EOL constraints (e.g., Mn computes 04 s). 

The system-level prediction results are shown in Fig. 3. Pre- 
dictions from the battery submodels are shown in Fig. 4. In 
this case, the motor temperatures reach a steady-state that is 
below the temperature threshold, so only the batteries have an 
impact on system EOL, which is the minimum of the EOLs 
predicted for the battery submodels. In particular, it is the 
first and fourth batteries (corresponding to Mn and .M 20 . 
respectively) that discharge the fastest, as shown explicitly in 
Fig. 4. The figures show the means of the predicted RUL 
distributions, the true RUL, RUL*, and an accuracy cone of 
a = 10% around it. In Fig. 3, we show both the system- 
level predictions using the distributed approach with Mn~ 
M 24 and the centralized approach using the global predic- 
tion model Mq- The global prediction model contains all 


the states, parameters, and constraints given in the previous 
subsection, minus the output constraints, and uses the com- 
manded wheel speeds (known a priori) as hypothesized in- 
puts. Since the currents are also known a priori, the system- 
level prediction can be decomposed, and the predictions made 
using the local submodels closely match those made using the 
global model, as shown in the figure. 

We use the relative accuracy (RA) metric (Saxena et al., 20 1 0) 
for prediction accuracy. Averaged over all predictions, RA is 
97.48% for the distributed approach and 98.74% for the cen- 
tralized approach. Using relative standard deviation (RSD) as 
a measure of spread, and averaged over all prediction points, 
RSD is 0.40% for the distributed approach and 0.43% for 
the centralized approach. The distributed approach is only 
slightly less accurate but has better precision. Here, both ap- 
proaches are very accurate since the system state-parameter 
estimates are very accurate, and there is only a small amount 
of error associated with assuming a constant average mo- 
tor current or wheel speed. Correspondingly, the prediction 
spread is relatively small because the uncertainty in the state- 
parameter estimate is very small. 

5.2.2. Friction Damage Progression 

We now consider a scenario in which for the front-left motor, 
there is nonlinear friction damage progression with v/fl = 
1 x 1CP 4 s. As a result of the continuously increasing fric- 
tion, the current drawn by the motor increases as well in order 
for the motor controller to maintain the same desired wheel 
speed. Hence, the total current drawn from the batteries is 
increased, and EOL occurs around 2 h. Because ipL is con- 
stantly changing, and in a way that is dependent on the motor 
state, it cannot be predicted a priori , and so cannot be used as 
a local input because the resulting predictions will not be ac- 
curate. Therefore, we require a submodel that estimates ipL , 
and we so employ submodels using as local inputs average 
values for the remaining motor currents, average commanded 
wheel speeds, and average rover translational velocity v and 
rotational velocity u. The prediction submodels for this case 
are shown in Table 3. For comparison, we demonstrate also 
prediction using M 17 -M 24 , and, for this strategy, at each 
prediction point the average value of current measured over 
the last minute is used as the future hypothesized value. Of 
course, this will not yield accurate results since future values 
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Table 3 . Prediction Submodels Using Commanded Wheel Speeds and Rover Velocities as Local Inputs 
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Figure 3 . System RUL prediction performance under nominal 
conditions with a = 0.1. 


Figure 4 . Individual battery submodel RUL prediction per- 
formance under nominal conditions with a = 0.1. 


of the current will actually be larger. Note that the predic- 
tion submodels used in this case do not correspond directly to 
those used for estimation. So, when constructing the estimate 
for M 25, for example, it takes the estimates from M\, Mg, 
and Mi3- 

The system-level prediction results are shown in Fig. 5 . Al- 
though the increased friction causes the temperature of the 
front-left motor to increase, it is still the batteries discharging 
that dominates the system-level EOL in this case. We show 
the predictions using M17-M 24 , M25-M32, and the global 
model Mo- For M25-M32 , average values of v = 0.5 m/s 
and uj = 0 rad/s are used. Here, the predictions using the lat- 
ter two approaches are virtually identical (the predictions us- 
ing M25-M32 are hidden under those for Mo), and fairly ac- 
curate. In contrast, as expected, the predictions using M\ 7 - 
M 24 are very inaccurate, and only converge towards the true 
RUL at the very end. This quite effectively demonstrates 
that, in this scenario, it is incorrect to use the front-left mo- 
tor current as a local input for predictions, since it cannot be 
predicted independently from the front-left motor submodel, 


and therefore a submodel that itself predicts this current is 
required to obtain accurate predictions. 

Here, RA averages to 58 . 95 % using M\ 7 -M2a, 94 . 24 % us- 
ing M25-M32, and 94 . 32 % using Mo- R SD averages to 
0 . 76 % using Mi 7 -M 2 4, 1 . 62 % using M25-M32 , 1 . 73 % 
using Mo. Here, we also observe an increase in prediction 
spread using the centralized approach with only a slight in- 
crease in accuracy over the distributed approach. 4 Overall, 
accuracy and precision are both decreased compared to the 
nominal scenario because there is more uncertainty in the 
state -parameter estimate, specifically, that dealing with the 
estimate of ujfl- This uncertainty in the state-parameter es- 
timate contributes to the additional uncertainty in the RUL 
predictions. 


4 The RSD for Mn~M2A is the lowest because those submodels do not 
include the motor friction component, so do not have the additional uncer- 
tainty associated with the estimation of the wear parameter. 
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Figure 5. System RUL prediction performance with friction 
damage progression with a = 0.1. 


6. Conclusions 

In this paper, we formulated the system-level prognostics 
problem and proposed a distributed solution based on struc- 
tural model decomposition. Using a four-wheeled rover as 
a simulation-based case study, we demonstrated the effec- 
tiveness of the approach. Most importantly, the distributed 
approach allows for, in many practical circumstances, the 
decomposition of the system-level prognostics problem into 
component-level prognostics problems and provides a mech- 
anism to merge local prognostics results into a system-level 
result. Further, since the local subproblems are independent, 
this allows component experts to focus on prognostics solu- 
tions for their components. However, we showed also that 
this approach is not always possible if accurate results are de- 
sired, since in some cases the prediction problem cannot be 
so easily decomposed, and it depends crucially on correct as- 
sumptions about what variables may serve as local inputs for 
the prediction problem. 

Although in this paper we focused on the model-based prog- 
nostics paradigm, our approach is flexible in that data-driven 
algorithms may be used also, once the local subproblems are 
defined. For example, in previous works, structural model 
decomposition was used to automatically design gray box di- 
agnosis models that were implemented using different data- 
driven techniques (for instance, state space neural networks 
in (Pulido, Zamarreno, Merino, & Bregon, 2012) or machine 
learning techniques in (Alonso-Gonzalez, Rodriguez, Prieto, 
& Pulido, 2008)). By decomposing the system-level prob- 
lem into independent subproblems through structural model 
decomposition, we can apply similar ideas to solve each prog- 


nostics subproblem by using the most appropriate technique, 
whether it is a model-based, data-driven, or hybrid approach. 

An important direction of future work is in algorithms for 
optimal placement of sensors for model decomposition, be- 
cause the level of model decomposition that can be achieved 
is dependent on the number of sensors and where they are 
placed. This results in the most efficient decomposition of 
the system-level prognostics problem. Current work also 
addresses combining the distributed prognostics framework 
with a distributed diagnostic approach for integrated diag- 
nostics and prognostics (Bregon, Daigle, & Roychoudhury, 
2012 ). 
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